Valley-dependent 2D transport in Si-MOSFETs 
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Motivated by interesting recent experimental results, we consider theoretically charged-impurity 
scattering-limited 2D electronic transport in (100), (110), and (lll)-Si inversion layers at low tem- 
peratures and carrier densities, where screening effects are important. We show conclusively that, 
given the same bare Coulomb disorder, the 2D mobility for a given system increases monotonically 
with increasing valley degeneracy. We also show that the temperature and the parallel magnetic field 
dependence of the 2D conductivity is strongly enhanced by increasing valley degeneracy. We ana- 
lytically consider the low temperature limit of 2D transport, particularly its theoretical dependence 
on valley degeneracy, comparing with our full numerical results and with the available experimental 
results. We make qualitative and quantitative predictions for the parallel magnetic field induced 2D 
magnetoresistance in recently fabricated high-mobility 6- valley Si(lll)-on-vacuum inversion layers. 
We also provide a theory for 2D transport in ultrahigh mobility Si(lll) structures recently fabri- 
cated in the laboratory, discussing the possibility of observing the fractional quantum Hall effect in 
such Si(lll) structures. 



I. INTRODUCTION 

Many semiconductor-based two dimensional (2D) elec- 
tron systems have an intrinsic valley degeneracy (g v ) in 
addition to the spin degeneracy (g s = 2). This valley 
degeneracy, which arises from the bulk band structure of 
the corresponding 3D material, is usually exact within 
the effective mass approximation, but is only approxi- 
mate in the experimental 2D systems where there could 
be small energy level splittings between different valleys 
(the so-called "valley splitting") 1 . If the valley splitting 
is "small" in some operational sense, the 2D system could 
be considered to have a total quantum degeneracy of 
g = g s g v with both spin and valley states being quantum 
degenerate. A well known example of valley degeneracy 
is graphene, which has a valley degeneracy of g v — 2 with 
equivalent Dirac cones at K and K' points of the Bril- 
louin zone^. Examples of semiconductor-based 2D elec- 
tron systems with valley degeneracy g v (> 1) are Si(100), 
(110), and (111) electron inversion layers in MOSFETs± 
as well as AlAs- and AlSb-based 2D electron systems^. 
Many other 2D systems (e.g. n-GaAs and p-GaAs 2D 
electron and hole systems) have no valley degeneracy 
(g v = 1)^. Valley degeneracy obviously has a profound 
effect on the electronic properties of the 2D system. The 
purpose of the current paper is a systematic theoretical 
investigation of the valley degeneracy effect on 2D elec- 
tronic transport properties, using Si-MOSFET structure 
based n-inversion layers as the specific system under con- 
sideration (since these typically have g v > 1) although 
our qualitative and analytical results would apply to all 
2D semiconductor systems (specifically Si-Ge 2D electron 
systems) with g v > 1. 

The valley-dependent electronic properties of a 2D (or 
3D for that matter) semiconductor system is best under- 
stood by considering the electronic density of states in the 
valley-degenerate ground state, which is given for the 2D 



(3D) system by: D(e) = grn/2nh 2 ( gv ^(2m) 3 / 2 /4:ir 2 h 3 ), 
where g = g s g v is the ground state degeneracy arising 
from both spin (g s ) and valley (g v ) degeneracies. The lin- 
ear proportionality of the electronic density of states with 
the valley degeneracy leads immediately to the following 
dependences of the Fermi wave vector (hp), the Thomas- 
Fermi screening wave vector (c[tf), and the Fermi energy 
(Ep) on the valley degeneracy in 2D and 3D systems: 
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Since D(e), kp, qxF, and ep all depend nontrivially on 
the valley degeneracy factor g v , all electronic properties, 
including 2D transport properties, depend nontrivially on 
the valley degeneracy. We note that the valley-dependent 
transport properties are in general nontrivial since g v en- 
ters independently through both qpp (and hence screen- 
ing) and kp (and hence scattering wave vector). Increas- 
ing (decreasing) the valley degeneracy enhances (sup- 
presses) screening through qTF, but at the same time 
it also affects Uf, increasing it with decreasing g v . 

Bulk Si has six equivalent conduction band minima 
located about 85% to the Brillouin zone boundary, thus 
making bulk Si a g v = 6 system. Each valley corresponds 
to an ellipsoid with anisotropic effective mass along and 
perpendicular to the symmetry axes. There are thus 
three possible Si-based 2D electron systems, depending 
on whether (100), (110), or (111) surface is used for cre- 
ating the 2D confinement. Within the effective mass ap- 
proximation these three distinct 2D Si systems have dif- 
ferent valley degeneracies: g v = 2 [Si(100)], 4 [Si(110)], 
and 6 [Si(lll)]. Uniaxial stress would lift the 4-fold or 
the 6-fold valley degeneracy of Si(110) or (111) system, 
making ground state of each a doubly degenerate g v = 2 
system similar to the Si(100) system (but, of course, with 
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a distinct effective mass for each 2D system). It is impor- 
tant to mention in this context the fact that essentially 
all Si-based 2D systems studied in the literature, inde- 
pendent of the surface orientation, have almost univer- 
sally always experimentally manifested g v — 2, even for 
Si(llO) and (111) systems, which nominally should have 
g v = 4 and 6, respectively. This is thought to be due to 
extensive random uniaxial stress universally present at 
the Si-Si0 2 interface in Si MOSFET structures, which 
pulls down two equivalent valleys compared with the 
other valleys, making both Si(llO) and Si(lll) MOS- 
FETs to have doubly degenerate ground states (g v = 2) 
similar to the Si(100) MOSFETs (but with different ef- 
fective masses) . An early experiment^ did manage to ob- 
serve a six-fold valley degeneracy in Si(lll)- SiC>2 MOS- 
FETs, but the system had very poor mobility and was not 
useful for the investigation of valley-dependent transport 
properties. Other than this one exception^, all Si-Si02 
2D MOSFETs invariably manifest g v — 2 ground state, 
independent of their surface orientation in sharp contrast 
to the effective mass approximation based expectation of 
9v = 6 (4) for Si(lll) [(110)] systems. 

An exciting new experimental development in the sub- 
ject, which is the direct motivation for our study, is the 
recent fabrication of very high quality Si(lll)-on-vacuum 
2DEG FET structures^, which exhibit g v — 6 ground 
state in agreement with the Si bulk band structure ef- 
fective mass approximation. Presumably the very high 
quality (without any interface strain) of the Si-vacuum 
interface, leading to very high mobility (~ 10 5 cm 2 /Vs), 
produce the expected 6-fold valley degeneracy. The fact 
that these Si-vacuum 2D FET systems also have very 
high mobility is consistent with the high quality of the Si 
surface leading to the g v = 6 Si(lll) 2D system. The ab- 
sence of a Si-Si02 interface may be the reason that these 
new Si-vacuum based 2D systems satisfy the expected 
g v = 6 effective mass approximation prediction. The 
absence of a real solid interface may simply enable the 
bulk effective mass approximation to be valid at the sur- 
face leading to the Si(lll) 6-fold degeneracy. Although 
similar 2D systems on Si(110)-vacuum system have not 
yet been made, it is reasonable to expect that the cor- 
responding Si(110) 2D system will have g v = 4 valley 
degeneracy. We note that these Si-vacuum 2D systems 
manifest, in addition to the expected valley degeneracy 
anticipated on the basis of the effective mass approxima- 
tion, also extremely high 2D mobilities because of the 
lack of random impurities in the oxide layer which ad- 
versely affects the mobility in the Si-Si02 2D systems. 

In particular, of course, the effective mass approxima- 
tion is not exact and breaks down at an interface 8 . Thus, 
even the double degeneracy (g v = 2) of the Si(100) 2D 
system [or of the Si(110) and (111) systems as observed 
experimentally] is only approximate and is lifted beyond 
the simple effective mass approximation leading to small 
(< lmeV) energy splitting between the two valleys. One 
can, therefore, think about an experimental single valley 
2D Si system where this ground state valley splitting is 



large enough so that the higher valley state is not occu- 
pied by electrons. Thus, in principle, the valley degen- 
eracy of Si-based 2D MOSFET or inversion layer system 
for any surface orientation can be thought to be a contin- 
uous variable ranging between 1 and 6 depending on the 
microscopic details of the interface. This is the approach 
we take in the current work where g v is assumed to be a 
free parameters to be determined experimentally. 

Motivated by the above considerations, we theoreti- 
cally consider valley dependent 2D transport in Si sys- 
tems assuming the valley degeneracy g v to be a free ra- 
tional variable — in reality, of course, g v can only be 
1, 2, 4, or 6 in Si 2D systems depending on the situa- 
tion. (In some cases, where the valleys are occupied by 
unequal number of electrons because of small valley split- 
tings, it may be useful to think of the valley degeneracy 
being a fractional number.) We address the density, the 
temperature, and the in-plane magnetic held dependence 
of 2D transport in the presence of a variable valley de- 
generacy. The applied in-plane magnetic field is parallel 
to the 2D system and is therefore assumed to only af- 
fect the spin degeneracy of the 2DEG since it gives rise 
to a Zeeman splitting between up/down spin levels. We 
find remarkably strong valley-dependence of 2D trans- 
port properties, and believe that the interesting physics 
of valley- dependent 2D transport should be investigated 
experimentally. The recently fabricated high-mobility 
Si- vacuum 2D electron systems&i should be particularly 
suitable in this context. 

We organize the rest of this article as follows: in sec- 
tion II we provide the detailed transport theory and a 
background giving a physical picture for why 2D carrier 
transport should depend strongly on the valley degener- 
acy, followed by the numerical results (III) . In section IV 
we provide our calculated mobility for the recently fab- 
ricated extreme high-mobility Si(lll) 2D samples. We 
conclude in section V with a discussion. 



II. TRANSPORT THEORY 

To calculate the density, temperature, and in- 
plane magnetic field dependence of 2D conductivity, 
cr(T, n, Bu), we use the Drude-Boltzmann semiclassical 
theory for 2D transport limited by screened charged im- 
purity scattering^. We assume that the 2D carrier con- 
ductivity is entirely limited by screened impurity scatter- 
ing, where the disorder arises from randomly distributed 
charged impurities in a 2D plane located at the interface 
between Si-Si02 (or Si- vacuum) and random background 
charged impurity centers (i.e. unintentional dopants) in 
the 2D layer itself. We neglect all phonon scattering 
effects as well as surface roughness scattering. In the 
low temperature limit (e.g. T < 10K regime of inter- 
est to us) phonon scattering is negligible for 2D elec- 
trons in MOSFET structures and the short-range surface 
roughness scattering at the Si-Si02 interface is only im- 
portant in the high-density limit (n > 10 12 cto -2 ) (the 
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roughness scattering can be neglected in Si- vacuum sys- 
tems). At low carrier densities (n < 10 12 cm -2 ) and 
at low temperatures (T < 10K) the 2D transport in 
Si-MOSFETs is dominated by the long-range Coulomb 
scattering by unintentional random charged impurities 
present at the Si-insulator interface and background un- 
intentional dopants inside Si. The background impurity 
density is low (~ 10 16 cto -3 ), but it may dominate all 
other scattering in high mobility Si- vacuum systems since 
the interface scattering is strongly suppressed due to the 
absence of the oxide layer. 

The 2D conductivity is given in the Boltzmann theory 

by 



J deD(e)^r(e) 



df(e) 
de 



(2) 



where D(e) = gm/ (2irh 2 ) is the density of states with 
total degeneracy g = g s g v and the carrier effective mass 
to, Vk = hk/m is the carrier velocity, e = (hk) 2 /2m is 
the usual parabolic 2D electron energy dispersion, /(e) 
is the Fermi distribution function, and r(e) is the energy 
dependent transport relaxation time. 

At T = we have a = ne 2 T(ep)/m, where n = gkp/An 
is the 2D carrier density with kp being the Fermi wave 
vector and ep the Fermi energy.. At finite temperatures 
we can express Eq. ([2]) by keeping the total carrier den- 
sity constant a = roe 2 (r) /to, where the energy averaged 
transport relaxation time (r) is given by 
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In the Born approximation the transport scattering time 
is given by considering screened charged impurity centers 
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where U z (q; z) is the bare Coulomb potential for electron- 
charged impurity interaction, q = k — k' is the momen- 
tum transfer, and Ni(z) is the random charged impurity 
density in the direction (z) normal to the 2D plane of 
confinement. The z = is the interface plane between Si 
and the insulator. In our model with two different kinds 
of impurity Ni{z) = Nf D + rii6(z - z ), where Nf D is 
the background 3D charged impurity density and ro^ is 
the 2D charged impurity density at zo from the inter- 
face. In Eq. (Ql e(q) = 1 — v(q)U(q) is the random phase 
approximation (RPA) 10 dielectric screening function due 
to the 2D electrons themselves, where v(q) is the 2D bare 
Coulomb interaction, and n(g) = n(g,T, By) is the 2D 
finite wave vector polarizability function depending on 
both temperature and spin polarization in the presence 
of finite parallel magnetic field. 

To understand a physical picture for why 2D carrier 
transport should depend strongly on the valley degener- 
acy we consider the conductivity for a strict 2D system 



at the zero temperature. However, in the realistic Si- 
MOSFET systems quasi-2D quantum form factor effects 
arising from the finite width of the 2D layer in z-direction 
must be included in V , '(q;z) and v(q). For this purpose 
the usual Howard- Fang variational function 11 is used in 
our numerical calculations. We will show that the quan- 
tum form factor effects are of considerable quantitative 
importance especially at low densities, where the system 
cannot really be thought of as an almost zero-width 2D 
layer. 

Using 2D RPA screening function at T — we have 
the scattering time for charged impurity centers 



_J_ = 2 ^A 2 g 2 /(g ), 
T{e F ) to g 

3/2 
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where go = gTf/2/cF oc g v (qrF — 9 /&b is a 2D 
Thomas-Fermi wave vector with effective Bohr radius 
as = h 2 K/me 2 where k is the background dielectric con- 
stant), and I(qo) is given by 



I{q ) w it + [2 + 41og(g /2)] qo for q 
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Thus in the strong screening limit (or qo ^> 1) we have 
T~ 1 (ep) = 7r2 ^^(f) 2 9v 2 > an d the conductivity be- 
comes a oc g 2 , i.e., the conductivity increases quadrati- 
cally with g v . In the opposite limit (i.e. weak screening 
limit, go < 1) we have r- l {e F ) = 2n 2 h r ^ (j) 2 q 2 , oc g v , 

and the conductivity becomes a oc g~ . In general, 
the effective screening (qo) of a 2D system becomes 
stronger as the density decreases. For Si(100) samples, 
qo fa A.7g 3 / 2 /y/h, where n = n/(10 10 cm~ 2 ). Thus, for 
the density regime of interest to us, n < 10 12 cm~ 2 , 
go > 1 and the transport of 2D Si-systems depends 
strongly (quadratically) on the valley degeneracy. 

At finite temperatures (T < Tp, where Tp = ep/kB is 
the Fermi temperature), the leading order correction to 
the conductivity is linear in temperature and given by- 



a(T) - cr 



1 - 



2go T 
l + qoT F 



(7) 



where <tq = a(T = 0). In the strong screening limit 
(go ^ 1) we- have 



Sa/ao oc —g v T/n, 



(8) 



where 5a = u(T) — do- On the other hand, in the weak 
screening limit we have 



8a I '(To oc 
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Thus for go 3> 1 the conductivity decreases linearly with 
valley degeneracy for fixed temperature and density, but 

for g <C 1 it decreases as g'J 2 ■ 

Since there are two different carrier components (spin 
up and down) in the presence of a finite parallel magnetic 
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field the total conductivity of the partially polarized sys- 
tem is given by a = er + + cr_, where a± = n±e 2 T±/m 
is the conductivity of spin up (+) and down (— ), re- 
spectively. n± is the carrier densities of spin state ±, 
and t± is the transport relaxation time of the spin 
up (down) state. To calculate the conductivity with 
screened charged impurities in the presence of parallel 
magnetic field, spin-polarization effects must be included 
in the polarizabilit y 12 ! 13 . When the parallel magnetic 
field is applied to a 2D electron system the polarizability 
becomes II tot (g) = H+(q) + H-(q), where n±(g) is the 
polarizability of the spin up (down) state and is given by 
at T = 



n±(g) - D F 



1 



^l-{2k$/q) 2 e{q - 2fc±) 



(10) 



where Dp = g u m/2'K 1 g v is the valley degeneracy factor, 
and kp are the Fermi wave vector of the spin up (down) 
state. 

For strictly 2D systems with zero thickness the spin 
polarization changes the screening function and this ef- 
fect gives rise to positive (negative) magnetoresistance 
in the strong (weak) screening limits 12,13 . In the strong 
screening limit (q ^> 1) a(B s )/a(0) » 1/4 and in the 
weak screening limit (qo -C 1) a(B s )/a(0) ~ 2, where 
cr(0) = <j{B = 0) and B s is the magnetic field for com- 
plete spin-polarization. Since the valley degeneracy af- 
fects the screening strength qo it is expected that the 
positive magnetoresistance is enhanced as g v increases. 
If we include the finite width confinement effect in the 
calculation the ratio becomes much smaller, especially at 
low density. In addition to the screening effects in the 
impurity potential, for the real systems with finite width 
confinement the orbital effects^ dominate over spin ef- 
fects at large magnetic fields. However we neglect the 
orbital effects in this paper, which has been considered 
elsewhere^. 



III. RESULTS 

Throughout this paper we use the following parame- 
ters: the dielectric constants of Si, SiC>2, and vacuum are 
Ksi = 11-7, Ksio 2 = 3-9, and n vac = 1, respectively, and 
the effective masses corresponding to the Si surface of 
(100), (110), and (111) arc m = 0.19m e , m = 0.28m e , 
and m = 0.30m e , respectively, where m e is the free- 
electron mass. In the absence of a parallel magnetic field 
the spin degeneracy g s = 2 is used everywhere. 

In Fig. [TJa) we show our calculated zero magnetic 
field conductivity of SiO 2 -Si(100) MOSFETs as a func- 
tion of valley degeneracy (g v ) for various electron den- 
sities n =1.0, 2.0, 5.0xlO n cm~ 2 (bottom to top) at 
T = 0, by assuming that the impurities with density 
rii = 3 x 10 n cm~ 2 are located at the interface (z = 0). 
As shown in the inset of Fig. Ufa) in logarithm scale, the 
calculated conductivity increases as cr{g v ) <x g™, where 
a increases as g v increases and approaches 2 which is 
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FIG. 1. Calculated conductivity for SiO 2 -Si(100) MOSFET 
system as a function of valley degeneracy (g v ) (a) for various 
electron densities n =1.0, 2.0, 5-OxlO 11 cm -2 (bottom to top) 
at T = and (b) for a density n =5. Ox 10 11 cm -2 and for 
various temperatures T = 0, 3, and 5K. Inset in (a) shows 
the same results of (a) in logarithm scale showing a ~ g% 
for large g v . Inset in (b) shows mobility at T = OK as a 
function of density for different valley degeneracies, g v — 1, 
2, 4, 6. The parameters corresponding to Si(100) are used 
and the impurity density of m = 3 x 10 11 cm -2 located at the 
interface (z = 0) is used. 



expected in the strong screening limit. In Fig. |Tfb) the 
conductivity is shown as a function of g v for a density 



n =5.0xlO u c 



and for various temperatures T = 0, 



3, and 5K with the same parameters as Fig. [lja). Wc 
see that the temperature dependence of conductivity be- 
comes stronger since the conductivity decreases linearly 
with g v for a fixed density and temperature in the strong 
screening limits [see. Eq. (|5|)]. 

In Figs. O and [3] we show the calculated resistivity 
= 1/cr) for Si0 2 -Si MOSFET systems assuming that 
all impurities are located at the interface z = 0. In 
Fig.[5]magnetoresistivity for (a), (b) Si(100) with g v = 1, 
g v = 2, respectively, (c) Si(110) with g v = 4, and (d) 
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FIG. 2. Magnetoresistivity of Si0 2 -Si MOSFET systems. 
Magnetoresistivity (p = 1/cr) of (a)(b) Si(100) (g v = 1, 
g v = 2), (c) Si(110) {g v = 4), and (d) Si(lll) (g v = 6) MOS- 
FET systems for various electron densities n =1.0, 2.0, 5.0, 
10.0 xlO 11 cm -2 (top to bottom). Impurities are located at 
the interface 2 = 0. Here p(0) is the resistivity at B = and 
B s is the magnetic field for complete spin-polarization, i.e. 
g a — 2 at B — and g 3 — 1 at B > B a . 
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FIG. 3. Resistivity of Si0 2 -Si MOSFET systems. Tem- 
perature dependent resistivity of (a)(b) Si(100) (g v = 1, 
g„ = 2) (c) Si(110) (g v = 4), and (d) Si(lll) (g v = 6) MOS- 
FET systems for various electron densities n =1.0, 2.0, 5.0, 
10. Ox 10 11 cm -2 (top to bottom). Here p(0) is the resistivity 
at T — and Tf is the Fermi temperature, Tf = EfI^-b- 
The same parameters as Fig. [2] are used. 



Si(lll) with g v — 6 MOSFET systems are shown for var- 
ious electron densities n =1.0, 2.0, 5.0, 10. Ox 10 n c77i -2 
(top to bottom). p(0) is the resistivity at B — and -B s 
is the magnetic field for complete spin-polarization, i.e. 
g s = 2 at B = and g s = 1 at B > B s . As expected the 
ratio p(B s )/p(0) is close to 4 for larger g s systems due to 
the enhancement of effective screening. Since the finite 
confinement effects the ratio p(B s )/p(0) is smaller than 
that of strict 2D system. Note that the results for Si(110) 
and Si(lll) with g v — 2 is very close to the results of 
Fig-H^b) even though we use the parameters correspond- 
ing to Si(110) and Si(lll). The results for Si(lll) with 
g v = 4 is almost identical with the results of Fig. [2jc) . 
The parameters corresponding to sample properties ex- 
cept the valley degeneracy have only a small effect on 
the ratio. The most important parameter determining 
the ratio is the screening strength go- For g v = 1 the 
effective screening becomes weak, and especially at high 
densities, qo < 1 is expected, as a consequence, the mag- 
netoresistance decreases as the magnetic field increases 
(see the result for g v = 1 and n = 10 12 cm~ 2 ). 

In Fig. [3] the temperature dependent resistivity of 
(a),(b) Si(100) (g v = 1, g v = 2, respectively), (c) Si(110) 
(g v = 4), and (d) Si(lll) (g v = 6) MOSFET systems 
are shown for various electron densities n =1.0, 2.0, 5.0, 
10. Ox 10 n cm~ 2 (top to bottom). Here p(0) is the re- 
sistivity at T = and Tp is the Fermi temperature, 



Tf = EfI^b- As shown in Eq. (jHJ at low tempera- 
tures T < Tp the metallic behavior (i.e., dp(T)/dT > 0) 
is strong for larger g v and at small densities. Thus, a 
stronger metallic behavior is expected for Si(lll) with 
g v = Qj£ For g v — 1 effective screening is weak, i.e., 
20 = ITF/Skp < 1, and strong-screening condition can 
only be satisfied at very low carrier densities. Thus, we 
expect rather weak temperature and field dependence of 
resistivity for g v — 1 except at very low densities. Fig. [3] 
shows the very interesting feature that the crossover tem- 
perature (i.e. the temperature at maximum resistivity) 
is almost independent of g v . Since Tp oc g v /n, the scaled 
crossover temperature (T /Tp) increases as the valley de- 
generacy increases for a fixed density, but the absolute 
crossover temperature is very close for all g s values. 

In Figs. |4] and [5] we show the calculated resistivity, 
p(B^) and p(T), respectively, for vacuum-Si MOSFET 
systems assuming that all impurities are located at the in- 
terface z — 0. In the calculation of resistivity for vacuum- 
Si MOSFETs we use the same parameters and the same 
impurity configuration for Si02-Si MOSFETS, except the 
insulating dielectric constant (Si02 versus vacuum, i.e., 
K Si0 2 versus K vac ), for direct comparison with the re- 
sults of Si0 2 -Si MOSFETs (Fig. and Fig. [3j. Since 
the screening strength qo is inversely proportional to the 
background dielectric constant at the same carrier den- 
sity both magnetic field dependent and temperature de- 
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FIG. 4. Vacuum-Si MOSFET systems. Magnetoresistivity of 
(a), (b) Si(100) (g v = 1, g v = 2), (c) Si(110) (g v = 4), and 
(d) Si(lll) (g v = 6) MOSFET systems for various electron 
densities n =1.0, 2.0, 5.0, lO.Ox IQ 11 cm~ 2 (top to bottom). 
All impurities are located at the interface 2 = 0. 



FIG. 5. Vacuum-Si MOSFET systems. Temperature depen- 
dent resistivity of (a)(b) Si(100) (g v = l,g v = 2), (c) Si(110) 
{g v = 4), and (d) Si(lll) (g v = 6) MOSFET systems for 
various electron densities n =1.0, 2.0, 5.0, 10. Ox 10 11 cm -2 
(top to bottom). Calculations are done with only interface 
impurities. We use the same parameters as Fig. [3] 



pendent resistivity for vacuum-Si MOSFETs is stronger 
than for Si0 2 -Si MOSFETs as shown in Figs. H and El 
However overall behaviors for both systems are very sim- 
ilar if we assume that the impurity configurations are 
identical. But in reality the impurity configurations are 
very different for these two systems. In the H-passivated 
Si-vacuum MOSFET it is expected that interface quality 
between vacuum and Si is much better with substantially 
less interface charged impurities than in Si-Si02- Espe- 
cially in high mobility vacuum-Si MOSFETs it is consid- 
ered that background unintentional 3D charged impurity 
is the most important scattering source. 

In Figs. [B] and [7] we show the calculated magnetic 
field and temperature dependent resistivity of vacuum- 
Si MOSFETs, p{B\\) and p(T), respectively, assuming 
that background charged impurities with impurity den- 
sity N? D = 2.3 x 10 16 cto -3 is the only scattering source 
(i.e. we set the interface impurity rii = 0). Even though 
the overall magnetic field and temperature dependences 
of resistivity are very similar to the results of Figs. 0] and 
[5] (which are calculated with interface charged impurity 
scattering) there are substantial differences between the 
two different impurity configurations. Comparing Fig. [5] 
with Fig. @] we find that the parallel magnetic field de- 
pendence of resistivity calculated with background 3D 
impurities is weaker than that with interface impurities 
(Fig. U). However the temperature dependence of resis- 
tivity calculated with 3D impurities is stronger than that 
with interface impurity (Fig. [5]) . Our predictions can be 



directly verified by carrying out parallel field measure- 
ments in Si- vacuum 2D systems. 



IV. HIGH-MOBILITY SI(lll) SYSTEM 

Very recent experimental work^ 6 . using ultra-clean H- 
passivated Si(lll)-vacuum 2D electron systems shows 
unprecedented high mobilities, approaching several hun- 
dred thousand err? jV s at low temperatures, correspond- 
ing to a momentum relaxation time (level broadening) of 
30 ps (~ 0.02 meV). These mobility numbers of these 
recent samples surpass the old Si(lll) mobilities^ by fac- 
tors of hundred and are comparable in quality (in terms of 
the momentum relaxation time and level broadening) to 
better-quality 2D GaAs electron samples where fractional 
quantum Hall phenomena typically manifest. Since the 
Si(lll) 2D system is a multi-valley system in contrast to 
the 2D GaAs system, the observation of the fractional 
quantum Hall effect in Si(lll) 2D system is a very inter- 
esting and potentially very important new development. 

In Fig [5] we show our numerical results (based on the 
theory given in section II) for the low-temperature mo- 
bility (defined simply as the conductivity divided by ne, 
(i = g /ne, where n is the carrier density and e is the elec- 
tron charge) of the Si(lll) 2D system, both for g v — 2 
and 6, as functions of temperature and density for sev- 
eral different impurity configurations. The density de- 
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FIG. 7. Vacuum-Si MOSFET systems. Temperature depen- 
dent resistivity of (a),(b) Si(100) (g v = 1, g v = 2), (c) Si(110) 
{g v = 4), and (d) Si(lll) {g v = 6) MOSFET systems for vari- 
ous electron densities n =1.0, 2.0, 5.0, 10. Ox 10 11 cm -2 (top to 
bottom). The same impurity configuration of Fig. [6] is used. 
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FIG. 8. Vacuum-Si(lll) MOSFET systems. The calculated 
density dependent mobility for different impurity configura- 
tions at T = 0.2K is shown (a) for g v = 2 and (b) for 
g v = 6. The black dashed lines represent the results with 



only 3D bulk impurities with Nf' 



2.3 x 10 1 



and 



other three lines represent results with only 2D impurities 
(no bulk impurities): z = OA with m = 2.0 x 10 10 cm -2 , 



z = 10A with m = 2.4 x 10 1 



and z — 100 A with 



rii — 21.0 x 10 cm~ . (c) and (d) show the calculated mo- 
bility with only 3D bulk charged impurities as a function of 
density for different temperatures T = 0.2, 0.5, 0.8, and 1.1 K 
(from top to bottom). In (c) g v = 2 and the bulk 3D impuri- 
ties of density Nf D = 4.5 x 10 15 cm -3 are used. In (d) g v = 6 
and the bulk 3D impurities of density iV 3D = 2.3 x 10 16 cm" 3 
are used. 



pendence agrees very well with the unpublished work of 
Kane 10 ', and shows that the disorder in this new batch of 
Si(lll)-vacuum samples is approaching the ~ 10 10 cm -2 
limit which is one to two orders of magnitude lower than 
the usual Si 2D samples^, thus explaining the very high 
mobilities of these new 2D systems. Our results also 
demonstrate that the strong screening by the g v = 6 
Si(lll) system, compared with the g v — 2 system, would 
lead to much stronger temperature dependence of the 
conductivity or the mobility, thus providing a clear means 
to distinguish the valley degeneracy. Our calculated den- 
sity dependence of the mobility is approximately consis- 
tent with the recent measurements^. 

Based on these calculation we can approximately 
estimate^ the expected activation energy of the 1 /3 frac- 
tional quantum Hall state of the Si(lll) 2D system (in 
these high-mobility samples) to be ~ 5 — 10K using 
the Zhang-Das Sarma models (and subtracting out the 
level broadening effect using our calculated mobility). 
The observation of the fractional quantum Hall effect in 
Si(lll) 2D system here is only possible because of the in- 
credibly high mobilities achieved through hydrogen passi- 
vation, and it is indeed a materials science breakthrough. 
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V. CONCLUSION 



The main purpose of this article is a systematic theo- 
retical investigation of the valley degeneracy effect on 2D 
electronic transport properties in Si-MOSFET systems. 
We calculate theoretically charged-impurity scattering- 
limited 2D electronic transport in Si(100), (110), and 
(111) inversion layers at low temperatures and carrier 
densities, where screened charged impurity scattering is 
important. The 2D mobility for a given system increases 
quadratically with increasing valley degeneracy, [i oc 
a/n oc gl, in the strong screening limit (go = ?T_F/2fci?) 
for the same impurity configuration. We also show that 
the temperature and the parallel magnetic field depen- 
dence of the 2D conductivity is strongly enhanced by in- 
creasing the valley degeneracy. All our results are valid 
only at carrier densities above which localization effects 
become important, but we estimate that the transition to 
the insulating state occurs well below 10 11 cm -2 density 
in the high mobility systems of our interesl^i^^ in the 
current work. 

We conclude by emphasizing our findings in both Si- 
Si0 2 and Si-vacuum MOSFETs. The parallel mag- 



netic field and the temperature dependence of the re- 
sistivity (at zero parallel field) manifest strong valley 
dependence regardless of impurity configurations. For 
g v = 1 both p(T) and p(B») show weak temperature 
and magnetic field dependence by virtue of weak screen- 
ing, and for g v = 6 due to strong screening (i.e. large 
90 = iTF/^kp) the resistivity shows both strong tem- 
perature and magnetic field dependence. Our finding 
of remarkably strong valley-dependence of 2D transport 
properties in Si-MOSFETs and the interesting physics 
of valley-dependent 2D transport should be investigated 
experimentally. Similar strong valley degeneracy depen- 
dence is also apparent in the many-body effects of 2D sys- 
tems, which have been studied elsewhere^. We have also 
provided detailed calculations for the valley- dependent 
transport properties of 2D Si(lll) systems as a function 
of temperature and density in very high-mobility low- 
disorder samples, commenting on the possible activation 
energy for the fractional quantum Hall effect in such ul- 
traclean Si system. 
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